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ABSTRACT 

We present a new method for directly determining accurate, self-consistent cluster 
lens mass and shear maps in the strong lensing regime from the magnification bias of 
background galaxies. The method relies upon pixellisation of the surface mass density 
distribution which allows us to write down a simple, solvable set of equations. We 
also show how pixellisation can be applied to methods of mass determination from 
measurements of shear and present a simplified method of application. The method is 
demonstrated with cluster models and applied to magnification data from the lensing 
cluster Abell 1689. 
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1 INTRODUCTION 

The possibility of reconstructing cluster lens mass distribu- 
tions from the magnification bias of background galaxies was 
first suggested by Broadhurst, Taylor & Peacock (1995) and 
first demonstrated by Taylor et al. (1998, T98 hereafter). 
They showed how a direct, local measure of the lens conver- 
gence, K = E/Ec, where E is the mass surface density and E^ 
is the critical surface density, could be obtained from knowl- 
edge of the lens magnification. In this way, one could mea- 
sure absolute surface mass densities, thereby breaking the 
"sheet-mass" degeneracy found in methods based on distor- 
tions of background galaxies (Tyson, Valdes & Wenk 1990, 
Kaiser & Squires 1993, Seitz & Schneider 1995). 

van Kampen (1998) and T98 have shown how one can 
extend magnification analysis into the strong lensing regime. 
By making reasonable assumptions about 7, the lens shear, 
they showed that one could place quite stringent bounds on 
K. In addition, T98 found an exact solution for the profile 
of axisymmetric lenses, although not for more general 2D 
cases. 

Inverse reconstruction methods based on maximum 
Hkelihood (Bartelmann et al. 1996) and majcimum entropy 
(Seitz, Schneider & Bartelmann 1998, Bridle et al. 1998) 
have gone some way in providing a unification of both shear 
and magnification information. Until now however, no direct 
method using only magnification has existed. 

In this letter, we show how to directly compute an accu- 
rate, self-consistent 2D distribution of k and 7 in the strong 
lensing regime from magnification. This direct approach has 
the advantage over indirect alternatives that uncertainties 
can easily be determined. The method is based on pixellisa- 



tion of the K distribution suggested by AbdelSalam, Saha & 
Williams (1998) who used it to estimate the mass of Abell 
370 from multiple images. We generalise the method further 
and also derive a simplified solution to the problem of esti- 
mating mass from shear, based on the approach of Kaiser & 
Squires (1993). 



2 RECONSTRUCTION OF k AND 7 

T98 showed how to estimate cluster surface mass using the 
magnification measured from the distortion in background 
galaxy number counts. Here our problem is to find an ac- 
curate method for reconstructing the surface mass density, 
given the magnification by an arbitrary lens. The inverse 
magnification factor at a given position in the lens plane is 



^-^ = 1(1-^)^ 



2| 
-7 I, 



(1) 



where k is the lens convergence and 7 is the shear. The shear 
can be decomposed into two orthogonal polarisation states, 
71 and 72, which are related to the lens convergence by 



71 



1 a-2/c,2 



3|)k, 72 = 9 ^8182^. 



(2) 



where di = d/dOi and is the 2D inverse Laplacian. The 
total shear is given by -y^ = 7i + 72 ■ One might expect that 
equation (0) could be solved iteratively by first estimating k, 
using this to calculate 7 and then updating the estimate of k 
using equation (Q) again. This proves to be highly unstable 
in the strong lensing regime however, rapidly diverging after 
only a few iterations (Seitz & Schneider 1995). 



© 0000 RAS 



L2 S. Dye & A. N. Taylor 



To find a stable solution to equation (^), we first pixel- 
lise the image. Following AbdelSalam et al. (1998), we can 
now write 



1,2 



(3) 



with summation implied over index m and where Km and 
7" are the pixellised convergence and shear distributions 
respectively. The transformation matrices, DY^" , are 



and 
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with the integration acting over the m* pixel, x = On — dm 
is the difference between pixels m and n which are assumed 
to be square in calculating these analytic expressions. Equa- 
tion (^) can now be written as the vector equation, 

1-2k + kGk -VA"^ ^0 (6) 

where is the A'^- dimensional vector of pixellised inverse 
magnification values, is the transpose of the vector k of 
pixellised convergence values and 1 is the vector (1,1,1,- -). 
The matrix G is the N x N x N matrix 



Gpqn — ^pn^qn 



(7) 



where Sij is the Kronecker delta, and summation is only 
over indices p and q. The parity of the measured inverse 
amplification A~^{0) is handled by V which fiips from being 
+ 1 outside regions bounded by critical lines to —1 within 
such regions. 

The amplification equation in the form of equation (^ 
is the first main result of this letter. We can now solve for k 
given a measured inverse amplification. Having solved for n, 
the corresponding shear distribution can then be calculated 
from equation (^). 



3 APPLICATION TO CLUSTER MODELS 

We apply the method to two types of idealised cluster mod- 
els. Starting with a predetermined cluster mass density dis- 
tribution, the corresponding shear distribution is derived us- 
ing Fourier methods (see for example, Bartelmann & Weiss 
1994). From these, the resulting magnification is calculated 
from equation (^) and then windowed to remove boundary 
effects. Using equation (^), we solve for k. 7 is then solved 
using equation (pf). A grid of 32 by 32 pixels is used in both 
models. 



3.1 Truncated Isothermal Sphere Model 

We first test the method with a simple truncated isothermal 
lens model. The pixellated mass distribution is laid down 
using K, oc (r-f ro)~^, where r is the radial distance from the 
centre of the sphere and ro is a constant. 



Figure ^ shows the k and 7 distribution from which 
the magnification distribution was calculated, the solved k 
and 7 distribution and the difference between them. The 
plotted distributions are smoothed from the underlying grid 
and the white dashes highlight the lens' critical line. The 
residuals are shown as percentage deviations from the true 
distribution. These are less than one percent for k over most 
of the grid which is negligible in comparison to the errors 
typically found in practice from background clustering, shot 
noise (T98) and the uncertainties resulting from use of local 
« estimators (see van Kampen 1998). The recovered shear 
distribution is more affected although still fares better than 
7 calculated from uncorrected Fourier techniques. The main 
contribution to these residuals is from boundary effects aris- 
ing from trying to recover a nonlocal shear in a finite area. 
Since much work has been carried out into the removal of 
such effects (see Squires & Kaiser 1996 and Seitz & Schneider 
1995 for example) which have little impact on the recovered 
we shall address the problem elsewhere. 

3.2 Dumb-bell Mass Model 

The method was also tested with a more general dumb-bell 
model. Magnification was determined in the same fashion 
as for the isothermal model, setting a negative parity inside 
the critical lines, shown by the white dashes in figure ^ 
Once again the residuals between the initial and solved k 
are typically less than one percent, while those for 7 are 
typically 10% and again come mainly from boundary effects. 



4 PRACTICAL CONSIDERATIONS 

We solve equation (^ with the hybrid Powell method (NAG 
routine C05PCF) . The number of equations needed to solve 
for K, is equal to the total number of grid pixels which can 
prove computationally intensive for especially fine grids. We 
find that this is not a problem for grid resolutions used to 
measure magnification bias in practice however. The 32 by 
32 grid of pixels used for the models in section ^ was solved 
in approximately one minute on an average workstation. The 
residuals exhibit no noticeable dependence on grid size. 

The Powell algorithm is an iterative process and there- 
fore requires an initial estimate of the solution to start from. 
The choice of the initial estimate turns out to be irrelevant. 
We have tried a wide range of initial distributions and even 
starting from a uniform distribution arrived at the same final 
solution. 

We have found that correct choice of pixel parity (espe- 
cially for low grid resolutions) is essential in order to achieve 
a sensible result. Inappropriate assignment of parities to pix- 
els manifest themselves, as one would expect, by k, being 
overestimated when a pixel is wrongly assumed to lie inside 
a critical line and underestimated in the reverse situation. 
This gives a means of checking whether critical line positions 
have been properly defined by looking for large discontinu- 
ities in the k distribution. Models with dual critical lines 
requiring dual parity flips have also been tested and we find 
that Hi can be recovered just as well. 

Finally, to ensure that the method does not break down 
with noisy data, we introduced a random noise term to the 
amplification. Errors in n resulting from noise in the inverse 
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Figure 1. Truncated Isothermal Sphere Model: The initial ft and 
7 used to form the magnification distribution from which the 
solved K and 7 are derived. Underlying grid dimensions are 32 by 
32. White dashes show the position of the critical line. Contours 
are linearly spaced and set at the same levels in both k plots 
and in both 7 plots. Residuals are expressed as percentages of 

(^init ^solved )/'^init 



amplification propagate as one would expect from equation 
(P). For an isothermal lens we recovered the expected result, 
5k = (5A/2A^, indicating that pixellisation does not lead to 
spurious noise properties. 



5 APPLICATION TO ABELL 1689 

We apply the method to the magnification data presented in 
T98 for the lensing cluster A1689. A 12 by 12 grid is used as 
the best compromise between shot noise in galaxy counts per 
bin and the resolution of the derived k map. Identification of 
the critical line was achieved by locating giant arc positions 
in the observed image. 

Figure shows the solved mass density and shear distri- 



Figure 2. Dumb-bell Model: Critical lines are shown as white 
dashes. Underlying grid dimensions are 32 by 32. Linearly spaced 
contours are set at the same levels for k and at the same levels 
for 7. 



bution. Comparison with the mass density map illustrated 
in T98 (figure 6) which was produced with the sheet k esti- 
mator shows very similar structure. We find that the value 
of K at the peak calculated here is approximately 10% lower 
than the peak value in T98 since the sheet estimator over- 
estimates K inside critical line regions. This has little ef- 
fect on the total integrated mass of A1689 found in T98. 
The 7 distribution is shown for completeness although un- 
doubtedly suffers from boundary effects typically found in 
the models. 



6 SHEAR ANALYSIS 

Having shown that pixellisation allows us to accurately re- 
construct surface mass densities from magnification data, 
we now apply it to shear analysis. Shear analysis exploits 
the idea that a given distribution of images of galaxies lying 
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Figure 3. A1689 solved convergence and shear distributions. 
Darker areas represent a higher distribution density. White dashes 
show the observed critical line. The plots are smoothed from a 12 
by 12 grid with north up and east to the left. 



behind a lensing cluster will, in the statistical mean, have 
regions of lens-induced correlations in image orientation and 
ellipticity. Measuring the quadrupole moments of individual 
galaxy images enables the construction of a map of the el- 
lipticity parameters, Cij (Valdes, Tyson & Jarvis 1983). The 
ellipticity parameters relate to the surface mass density and 
shear via (Kaiser 1995, K95 hereafter): 



1 



71 
72 



72 

-71 



(8) 



One way of solving this for k in the weak lensing regime is 
to follow the approach of Kaiser & Squires (1993). Gener- 
alisations of this to the strong regime have been made by 
K95. One would have hoped that an alternative to such ap- 
proaches would be to pixellise equation (^ and use equation 
(^) to solve it by matrix inversion. However, the resulting 
matrix equation is ill-conditioned, since the matrix D™" is 
singular and _D™" is itself ill-conditioned. Instead, we show 
a new, simplified expression for the solution to Kaisers' el- 
lipticity equation and then pixellise it. 
Starting with the equation (K95), 



and using equation one can show that 
di ln(l ~ k) = ~dj \n{5ij + eij). 



(9) 



(10) 



The term on the right hand side is obtained from the defi- 
nition. 



ln(/ + B) = B + + + ■ 



(11) 



where / is the identity matrix and B is an arbitrary square 
matrix. Using this expansion and collecting even and odd 
terms we find. 



--ln(l-e2)5., + 2 



2 Vl-e^ 



(12) 



where = e\-\- e\ and Ci — 7^/(1 ~ k). This result requires 
that e < 1. Inserting equation (^) into the magnification 
equation (jl|) we find 

A'' ^\{l~K)\l~e% (13) 

Hence the parity changes when e > 1. Since eij and e^^ 
are observationally indistinguishable and fiip from one to 
another whenever there is a parity change, we can satisfy 
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Figure 4. Reconstruction of k from the ellipticity parameters. 
Contours are at the same levels in both n plots. The distortion 
field is illustrated by plotting the apparent shape of an intrinsi- 
cally circular background object. 



the criterion e < 1 just by noting the critical line positions 
and inverting the ellipticity matrix when one is crossed. 

Finally, inserting equation (jl^) into equation (|lo[), and 
solving for k we find the pixellised solution is 



o2^l/2 



exp 



where 



■In 



1,2 



(14) 



(15) 



Equation ( p^ ) is the second main result of this letter. 
We can directly calculate k. given a measured ellipticity field. 
Figure B shows the results of reconstructing k using equa- 
tion ( p,^ for the dumb-bell model. The ellipticity parameters 
are calculated from equation (^) using the k and 7 distribu- 
tion. We normalise the reconstructed k to both peaks in the 
initial k distribution. The residuals, again being dominated 
by boundary effects, show that reconstruction is possible to 
within approximately 10 percent across the field of view. 



7 SUMMARY 

We have outlined a method for directly calculating accurate, 
self-consistent surface mass density and shear distributions 
from the lens amplification and critical line positions. The 
method has been demonstrated with the isothermal sphere 
and dumb-bell cluster models. We find it reconstructs the 
surface density to within a percent over most of the field 
of view. The reconstruction of the shear pattern only has 
fractional accuracy of a few tenths due to boundary effects. 
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We have applied the method to magnification data from 
Abell 1689, and reconstructed its surface mass and shear 
distribution. 

We have also found a simplified solution to the prob- 
lem of estimating surface mass density from galaxy elliptici- 
ties. This approach puts the calculation of surface mass from 
shear and magnification on an equal footing, and we shall 
investigate the combined analysis elsewhere. 
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